Exploring Quantum Phase Transitions with a Novel Sublattice Entanglement Scenario 
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We introduce a new measure called reduced entropy of sublattice to quantify entanglement in spin, 
electron and boson systems. By analyzing this quantity, we reveal an intriguing connection between 
quantum entanglement and quantum phase transitions in various strongly correlated systems: the 
local extremes of reduced entropy and its first derivative as functions of the coupling constant 
coincide respectively with the first and second order transition points. Exact numerical studies 
merely for small lattices reproduce several well-known results, demonstrating that our scenario is 
quite promising for exploring quantum phase transitions. 
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Quantum entanglement, a key concept in quantum in- 
formation theory that has no classical counterpart 0, 
has been recognized to play also an important role in the 
study of quantum many-particle physics Q. However, 
how to quantify appropriately entanglement in many- 
body systems has been a challenging question for a long 
time. So far, a number of interesting theoretical at- 
tempts have been made on various notions of entangle- 
ment in several one-dimensional (1D) correlated systems 
as well as their possible links to quantum phase tran- 
sitions (QPTs) that oceur at absolute zero temperature 
with the change of coupling parameters |ÏÍ 1110 SES Oi- 
Most of these studies focused mainly on the spin chains, 
where the spin-spin concurrence pj was adopted to de- 
scribe (two-particle) entanglement, which can be calcu- 
lated from correlation functions Since the correla- 
tion function decays rapidly, the concurrence is nonzero 
only between two closer sites and thus the information 
on QPTs abstracted from the concurrence is rather lim- 
ited 0] ■ Scaling behavior of entanglement between a 
block of contiguous spins and the rest of the system in 
a spin chain was also studied around the quantum criti- 
cal point (3 , but the idea of the contiguous-spins entan- 
glement is merely applicable to few 1D models without 
classifying the transition order. Therefore, it is of sig- 
nificance to attempt other more applicable measures of 
entanglement for correlated systems so that the nonlocal 
nature of QPTs can be better captured. 

In this paper, by introducing a well-defined reduced en- 
tropy of sublattice as a measure of entanglement, which 
is an extensive quantity, we have found a distinct connec- 
tion between quantum entanglement and QPTs in such 
an intriguing way: the local extremes of reduced entropy 
and its first derivative as functions of the coupling con- 
stant coincide respectively with the first and second order 
transition points of QPTs. This finding is remar kable as 
(i) it implies that the reduced entropy of sublattice plays 
a crucial role in QPTs, similar to that of the thermal 
entropy in classical phase transitions; and (ii) the analy- 
sis of reduced entropy enables us to identify and classify 
unambiguously QPT points for a wide class of systems. 

While we have studied more than ten strongly corre- 
lated systems with the same conclusion, here we first an- 



alyze the best understood transverse-field quantum Ising 
chain as a benchmark, and then present our main re- 
sults for three types of two-dimensional (2D) spin-1/2 
systems and one type of electron system: a coupled 
dimer antiferromagnet (Fig.la), a frustrated J\-Ji model 
(Fig.lb), a checkerboard antiferromagnet (Fig.lc), and a 
1D half-filled extended Hubbard model (Fig.ld). The 
corresponding Hamiltonians read 



H a — 2_j ^ij Si Sj i 



(1) 



(ij)a i (ij) 

where Si is the spin operator and a is a model label. The 
coupling parameters in all models are explicitly shown 
in Fig.l. For the coupled dimer antiferromagnet model 
(a = a), the anisotropy parameter is denoted by A. 
J\ and Ji correspond to the nearest-neighbor and next- 
nearest-neighbor exchange couplings in the 2D frustrated 
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FIG. 1: Schematic configurations of the four studied strongly 
correlated systems. a) Coupled dimer antiferromagnet, b) 2D 
frustrated J\-Ji system, c) checkerboard antiferromagnet and 
d) 1D extended Hubbard model. Two subsystems are denoted 
as X (crosses) and • (sòlid symbols) in the figures. 
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Heisenberg Hamiltonian (a = b). In the checkerboard 
antifcrromagnet spin system (a — c), the Hamiltonian 
contains nearest-neighbor couplings along the nearest- 
neighbor coupling J and the diagonal link J x . As for the 
1D half-filled extended Hubbard model, Cj a and Cj a arc 
creation and annihilation operators of electron with spin 
o at site j, respectively, U and V represent the on-site re- 
pulsion/attraction and nearest-neighboring interaction. 

In an earlier paper, Vidal et al. ,7j studied the ground 
state entanglement between a block of contiguous spins 
and the rest of the system. Their approach is, however, 
not applicable to the above four systems and many oth- 
ers. In the present work, it is in fact crucial to choose 
one size-N sublattice Rn reduced from the whole lattice 
in such a preferable way (i) the number of connecting 
bonds between the chosen sublattice and the rest (-B/v) 
of the system is maximized so that the correlation effect 
between the two parts can be best revealed; and (ii) the 
structure of sublattice Rn should be the same as (or, if 
impossible, most similar to) that of the original lattice, 
in the same spirit of the real-space renormalization idea. 
Under such choice, a new appropriate measure of entan- 
glement for correlated systems, i.e. the reduced entropy 
of sublattice, is naturally given by Sn — — tr (/?n l°g 2 Pn), 
where pn is the reduced density matrix by the partial 
trace of the density matrix over Bn Q . Being different 
from some conventional quantities (e.g., the spin suscep- 
tibility) which arc only related to a two-body correla- 
tion function, this reduced entropy of sublattice contains 
all the contributions including two-body, three-body, and 
even iV/2-body correlation functions. Since this reduced 
entropy can capture essentially the nonlocal nature of 
many-particle correlation, particularly around the tran- 
sition points of QPTs, we naturally anticipate that it 
should have a close connection to QPTs, as we will elab- 
orate later. In the numerical calculations (except for the 
first example addressed below), we employed Lanczos al- 
gorithms with periòdic boundary conditions to calculate 
exactly the ground state ^g), from which we obtained 
Pn = tre N |í'g)(í'g| and then Sn, with the lattice size 
N and choice of the sublattice Rn (the lattice of X- 
points (crosses)) used in our calculation being illustrated 
in Fig.l. 

As a benchmark illustration, we first analyze an ex- 
actly solvable model intensively studied before 0, 0, , 
i.e., the transverse-field quantum Ising chain with the 
Hamiltonian as H = -£Í=i (SfSf +1 + XSf). It is well 
known that this system undergoes a second order phase 
transition between an ordered ferromagnetic state and a 
quantum paramagnetic state at A c = 1. As depicted in 
Fig. 2 (with the even sites sublattice), Sn monotoni- 
cally decreases as A changes from to 2, while the first 
derivative d(SN/N)/d\ shows explicitly a local minimum 
at A = A c , indicating clearly that this local extreme cor- 
responds to a critical point of the second-order QPT. 
This one-to-one correspondence between the reduced en- 
tropy and QPTs has also been tested and justified in 
the study of the Anderson localization model describ- 
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FIG. 2: The ground state renormalized reduced entropy 
Sn/N (a) and its first derivative versus the external param- 
eter A (b). The curves correspond to dififerent lattice sizes 
N= 10,40,160,1280. 



ing disorder-driven metal-insulator transition In 
fact, from our numerical results including those to be 
addressed later/elsewhere, we state more firmly that the 
extreme of the first derivative of reduced entropy in a fi- 
nite system corresponds to the transition point of second- 
order QPTs. Moreover, as shown clearly in Fig. 2, as 
the lattice size increases from N = 10 to 1280, the de- 
termined critical point is unchanged, and the curves of 
the renormalized first derivative (and reduced entropy) 
converge very rapidly (for N > 160, all curves merge to 
one within the present resolution) |l2| . Therefore, our 
reduced entropy of sublattice indeed captures essentially 
the nonlocal naturc of many-particle correlation and re- 
veals QPTs with negligible (or very weak) finite size ef- 
fect, which is a remarkable advantage in exploring QPTs 
in other complicated systems. 

We now use the above scenario to look into a coupled 
Heisenberg dimer Hamiltonian with two spins per unit 
celi, a typical model for illustration of novel QPTs in the 
literature |Ï3j. The transition may be tuned by varying 
a dimensionless parameter A (0 < A < 1). As shown in 
Fig.l (a), the bonds represented by the sòlid- line form 
the coupled dimers while the dashed-line represents the 
coupling between the neighboring dimers. In the case of 
A< 1, the ground state is a paramagnetic state consist- 
ing of a product of decoupled dimers, in each of which 
the spins pair into a valence bond singlet. At A=l, the 
system corresponds to the square lattice antiferromag- 
net; its ground state has a long-range magnètic Neel or- 
der. Obviously, these two limits cannot be connected 
continuously and a transition point is expected between 
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FIG. 3: The ground state renormalized reducec 
Sn /N and its first derivative versus the external p: 
a) Coupled Dimer Antiferromagnet, b) 2D frustra 
system. 



them. Recently, a field theory study indicated uw** 
transition between the two phases is of second-order 0] , 
and the critical point was estimated around 0.5 from a 
Monte Cario simulation [14| . For the sublattice specified 
in Fig.l(a), the reduced entropy Sn and its first deriva- 
tive as functions of coupling A are shown in Fig. 3(a). 
Sn monotonically decreases as A changes from to 1, 
while the first derivative d(SN/N)/d\ shows a clear lo- 
cal minimum around A ~ 0.4, which is now identificd 
to be a critical point of the second-order QPT, in good 
agreement with the well-accepted result. 

For the 2D frustrated J\-Ji Heisenberg model 
(Fig.l(b)), it is known that the ground state of the 
model possesses the conventional Neel long range or- 
der at J2/J1 -C 1 and has a collinear antiferromagnetic 
quasi- long range order if J 2 /Ji ~ 1. The reduced en- 
tropy and its first derivative are illustrated as functions 
of J2/J1 in Fig. 3 (b). By analyzing the first derivative 
cI(Sn /N) I d{Ji I Ji) as a function of J2/J1, we locate two 
local extreme points at J%/J\ ~ 0.37 and 0.62, respec- 
tively. According to the above criterion, we identify two 
critical points in the system. Besides, there is also a lo- 
cal maximum of Sn at J2/ Ji=0.5, which is predicted 
to be another first-order QPT point according to our 
earlier study |ïfij and many new results including those 
to be presented later. Interestingly, there is indeed a 
general consensus in previous studies that three distinct 
QPT points are present around J%j J\ ~ 0.38, 0.5 and 
0.6 [TE El EHH- In view of the small 4x4 cluster 
used in our calculation, the good agreement of our re- 
sults with the above consensus is very impressive. 

Next we consider one of the most frustrated 2D anti- 
ferromagnets, the checkerboard antiferromagnet, which 
was investigatedby some analytical and numerical ap- 
proaches before [ïjj . Its ground state has Neel long-range 
order at J x /J <C 1, whereas it corresponds decoupled 
Heisenberg chains when J x /J> 1. At J x /J = 1, it is 
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FIG. 4: (a) The ground state renormalized reduced entropy 
Sn/N (sòlid line) and its first derivative (dotted line) versus 
the external parameter J x / J for the 2D checkerboard anti- 
ferromagnet model, (b) The phase boundaries for the 1D 
extended repulsive Hubbard model. 



a valence bond crystal in singlet plaquettes. These three 
phases are separated by two quantum critical points at 
(J x /J)c ~ 0.80 and 1.10, respectively. A recent field 
theory study reported a new candidate phase diagram 
with additional intermediate ordered phase [2(|. In Fig. 
4(a), the calculated reduced entropy and its derivative 
are plotted as functions of J x / J ■ As we can see, there 
exist a local maximum in Sn and two local extremes in 
its first derivative located respectively at J x /J ~ 0.93, 
J x / J ~ 0.82, and 1.02. Based on our scenario, we predict 
that there exist eventually four different phases divided 
by three transition points. The same technique is also 
applied to study the hard-core bosons on the checker- 
board lattice, revealing its possible superfluid to super- 
solid transition |2l| . 

A typical electron example is the 1D extended Hub- 
bard model at half filling (Fig.l(d)). It is well known 
that there are two classical phases: charge-density- 
wave (CDW) (U/V <C 1) and spin-density-wave (SDW) 
(U/V 3> 1) for positive U. The classical first-order 
phase boundary lies at U/V—2 22]. By employing our 
approach, we can identify accurately this classical first- 
order boundary from the local màxima of the reduced en- 
tropy. On the other hand, a recent theoretical study pre- 
sumably indicated the presence of a new quantum phase 
bond-order-wave (BOW) at the verge of the line U /V=2, 
with the QPTs between SDW/BOW and CDW/BOW 
appearing to be the second-order (2^. By analyzing the 
local extreme of the first derivative d(SN/N)/dV at fixed 
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U, two critical points close to the line U/V=2 show up. 
Thus we incline to predict that two new quantum phases 
emerge around the classical phase boundary where the 
residual CDW and SDW correlations may still survive in 
the presence of quantum fluctuations, with the two new 
phase boundaries being the second-order (Fig.4(b)). By 
increasing the lattice size (from N=6 to 10), the shifts 
of local extreme lines of d(SN/N)/dV towards the line 
U/V=2 can be clearly observed (not shown here). Thus 
the phase region for BOW becomes narrower, which is 
consistent with a previous investigation 23]. It is no- 
table that the well-known second-order superconducting 
phase transition is also revealed with the same analysis 
for a negative U. 

For the above numerical finding of the relationship be- 
tween the reduced entropy of sublattice and the tran- 
sition points, particularly from the first example, we 
may heuristically undcrstand its origin: as N — > oo, the 
shape of the peaks or dips around the transition points 
(in Figs. 3 and 4) would become sharper |24|, namely 
d 2 (S N /N)/d 2 X (or d(S N /N)/dX) would be discontinuous 
at the second (first) order transition points, resembling 



to the cases in thermodynamic phase transitions. 

Many QPTs are yet to be clearly understood due 
to the competing nature of highly degenerate quantum 
statcs 9J. In the sense that the reduced entropy for a 
small system may capture essentially the nonlocal fea- 
ture of entanglement in a real system, this appropriate 
measure of entanglement enables us to establish its clear 
connection to transition points of QPTs. Therefore, com- 
paring with the conventional wisdoms, it is not compul- 
sory for us to deal with a large system to make the scal- 
ing analysis. In our scenario, exact numerical results for 
a quite small correlated system are able to disclose rele- 
vant information about QPTs. This is indeed a unique 
and superior advantage of our approach, making it be 
quite promising for the future exploration of QPTs in 
various physical systems, particularly if they fail to be 
charactcrized by the conventional approach. 
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